Substrate influence on the plasmonic response of clusters of spherical nanoparticles 
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The plasmonic response of nanoparticles is exploited in many subfields of science and engineering 
to enhance optical signals associated with probes of nanoscale and subnanoscale entities. We develop 
a numerical algorithm based on previous theoretical work that addresses the influence of a substrate 
on the plasmonic response of collections of nanoparticles of spherical shape. Our method is a real 
space approach within the quasi-static limit that can be applied to a wide range of structures. We 
illustrate the role of the substrate through numerical calculations that explore single nanospheres 
and nanosphere dimers fabricated from either a Drude model metal or from silver on dielectric 
substrates, and from dielectric spheres on silver substrates. 



I. INTRODUCTION 

Currently there is great interest in the use of the plas- 
monic response of tailored metallic substrates and other 
structures for the purpose of enhancing electric fields of 
laser beams in their near vicinity. Enhancements with 
origin in the excitation of collective plasmon modes can 
increase the field intensity by many orders of magnitude 
in the near vicinity of diverse systems. This phenomenon 
was first explored in the context of surface enhanced Ra- 
man scattering (SERS), wherein it was found that the 
Raman cross section of pyridine adsorbed on electro- 
chemically roughened Ag surfaces can be enhanced by 
approximately six orders of magnitude relative to that 
realized for pyridine in solution^-. The field has evolved to 
the point where the Raman spectrum of single molecules 
can be detected through use of plasmon enhanced Raman 
probeJ2E, Plasmonic enhancements can be used not only 
in the context of Raman spectroscopy, but more generally 
to enhan ce t he cross section of diverse nonlinear optical 
processe^^. 

In the theoretical literature, one finds numerous stud- 
ies of the pla smo nic response of isolated nanoparticles of 
diverse shap^^ along with metallic arrays of nanopar- 
ticleiPal. So far as we know, virtually all such dis- 
cussions explore nanoparticles and their arrays in free 
space^" Treatments of the free space response are ap- 
propriate for clusters of nanoparticles in solution, but 
commonly one is interested in particles and particle ar- 
rays on substrates. Then an issue is the influence of the 
nanoparticle-substrate interaction on the plasmonic re- 
sponse of the nanoparticles that reside on it. Papers ad- 
dressing particle-substrate interactions include the work 
of Yamaguchi et al.'^^, which discussed particles above 
substrates in the dipole appr oxima tion. Work done by 
Ruppirfi^ and by Noguez et al! ^^ * ^^ ! also deal with sphere- 
substrate interactions, but only for a single nanoparticle. 
Mayergoyz et al. have studied the plasmon eigenfrequen- 
cies of nanosphere dimers and also cylindrical structures 
on a substrat^. Moreover, a recent study on the plas- 



monic response of cubical nanoparticle dimerJi^ reports 
on the dimer-substrate interactions in the SERS context. 

Since the early 1970s, Bedeaux and Vlieger et al. have 
conducted numerous theoretical and numerical studies 
on the effects of particle-substrate and particle-particle 
interactionJi^ These studies have been concentrated 
around spherical or spheroidal particles on top of a sub- 
strate, or truncated particles of such shapes on a sub- 
strate (used to model a finite contact angle) . In this work, 
the particle-substrate interactions were taken into ac- 
count to high multipolar order, while the particle-particle 
interactions were only calculated to dipolar or quadrupo- 
lar order, since their main c oncern was systems of low or 
finite particle coverag^^^E^. 

More recently, numerical studies based on Bedeaux and 
Vlieger's work have been carried out by Simonsen, Laz- 
zari and co-workers for the purpose of in-situ inversion of 
experimental optical spectra obtained from growing thin 
granular metal films^^. 

In this paper, we present a description of the influence 
of a substrate on the plasmonic response of non-periodic 
nanosphere arrays; through use of the Bloch theorem one 
may address periodic systems as well. We employ the 
quasi-static description of the response of the system. 
This proves adequate for objects whose linear dimensions 
are small compared to the wavelength of Iigh1p2l. In con- 
trast to previous work, we consistently take into account 
higher order interactions between the nanospheres. 

After we describe the formalism, we turn our attention 
to calculations that explore the influence of the substrate 
on the response of nanospheres and nanosphere dimers. 
Of interest is the discussion of "hot-spot" -regions where 
at selected excitation frequencies one realizes very large 
field enhancements by virtue of the excitation of collec- 
tive plasmon modes. For the case of two spheres in free 
space that are nearly in contact, one realizes a ho t-spot 
at the point of closest contact between the sphere^^aiSH 
In this paper, for a nanosphere dimer near a dielectric 
substrate we find "moving hot-spots" . A small change in 
excitation frequency can cause the hot-spot to move from 
the point of nearest contact between the spheres, to the 
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"south poles" of the spheres — the points on the spheres 
closest to the substrate. In recent work, two of the 
authors have discussed moving hot-spots in nanosphere 
clusterj2ll. 

The present study illustrates the role of the substrate 
in creating new hot-spots. We find that if a dielec- 
tric sphere is in close proximity to a "plasmonic active" 
metallic substrate, the region around the south pole of 
the dielectric sphere becomes a hot-spot. A spatially 
localized "potential well" that can trap substrate plas- 
mons is formed just under the dielectric sphere. Also, if 
a metallic sphere is placed close to a dielectric substrate, 
we find a collective plasmon localized near the south pole 
of the sphere. Thus, the interaction of nanospheres and 
structured arrays of such objects placed on substrates 
creates new hot-spots that can be exploited in diverse 
non-linear optical spectral probes of nanoscale and sub- 
nanoscale matter. 

In this paper, Sec.[ll]is devoted to setting up a formal- 
ism that may be applied to any non-periodic structure 
of spherical nanoparticles that are located on, or near, 
a substrate, and Sec. [Ill] presents the results of our nu- 
merical studies of isolated nanospheres and nanosphere 
dimers on substrates. Section [Tv| contains concluding re- 
marks. 



II. THEORY 

Even if the numerical calculations to be performed in 
this paper will focus on one or two nanoparticles, we 
will however present a more general formalism valid for 
a cluster of N nanoparticles. For the case of the dimer, 
the geometry is illustrated in Fig. [T] The substrate is 
located in the half space z < and it is characterized by 
the dielectric function £_(a;). The region above the sub- 
strate, z > 0, is assumed to be a non-absorbing dielectric 
characterized by the dielectric function e+{uj). 

We consider a system consisting of N non-overlapping 
nanospheres, located at arbitrariy positions. For each 
such sphere we embed a coordinate system Sj, j = 
1,2, ... N, so that the origin of Sj is located at the center 
of sphere j. With each coordinate system Sj, we asso- 
ciate a position vector Vj = {rj,9j,4>j). 

Our interest is in nanosphere arrays whose extent is 
small compared to the wavelength of light, so the electro- 
static approximation suffices to describe the electric fields 
in its vicinity. Under this assumption, the Maxwell's 
equations are equivalent to the Laplace equation. Thus 
our task is to solve Laplace's equation for the electro- 
static potential ip, 

= 0, 

subject to the appropriate boundary conditions on the 
surface of each sphere, and at the interface between the 
substrate and the rest of the system. As usual, the elec- 
tric field is given bjJ^l E = —Vij}. 
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FIG. 1. (Color online) An illustration of the system we con- 
sider in this paper, for the case where we have a nanosphere 
dimer. The substrate occupies the half space z < 0, /i is the 
distance between the "south poles" of the spheres and the 
substrate, and d is their surface-surface separation. Sphere 
j of the dimer has dielectric function Ejiuj) and radius aj. 
The medium above the substrate has dielectric function e+(t<j) 
while that of the substrate is £-{uj). The two black dots rep- 
resent schematically image multipoles in the substrate seen 
by an observer in the half space 2 > 0. 



We will assume that a spatially uniform electric field 
of angular frequency lo is applied to the system and 
we analyze its response to this field. In what follows all 
dielectric functions that enter the analysis are the com- 
plex dielectric functions appropriate for the frequency lo, 
though we suppress any explicit reference to in what 
follows. Hence, the electrostatic potential in the half 
space z > can be written as 
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where is the electrostatic potential produced by the 
polarization charges in sphere j and V'j the potential pro- 
duced by its image, located in the half space z < 0. In 
the substrate (z < 0), the electrostatic potential takes 
the form 
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where ipj is the electrostatic potential of sphere j as seen 
by an observer in the region z < 0, and Eq is the applied 
field in the substrate. The various single sphere potential 
functions that enter Eqs. ([I]) and ^ may be expanded in 
the spherical harmonics. Using the shorthand notation 

/,„ = L;=o Lm=-i we have 



I 



(3a) 
(3b) 
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and 

V'J(^.) = E 4t^''^7'^'Yrio,,<Pj), (3c) 

Im 

where the various Aim and Bim are expansion coefhcients 
to be determined, and aj refers to the radius of sphere 
j. The symbol Fj™ refers to the spherical harmonic func- 
tions as described by Ref. [55] As discussed in Refs. [T7] 
and im the coefficents a\^^ and a\^^^ are related to 

^hn through the boundary conditions at the interface 
z = 0. Simple image arguments supply the relatio n Ije- 
tween these quantities. In particular, one finds thalP^E^l 



Ref. ini to incorporate interactions of periodic structures 
with a substrate. One then has plasmon normal modes 
characterized by a wave vector fcy parallel to the surface; 
one encounters only {L + 1)^ coefficients in this case, be- 
cause the expansion coefficients of different nanoparticles 
are linked by the Bloch theorem. The quasi-static limit 
developed in this paper can be applied to the descrip- 
tion of collective excitations whose wave vector is large 
compared to w/c, with cj the angular frequency of an ex- 
citation of interest, and c the velocity of light in vacuum. 



III. RESULTS AND DISCUSSION 
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Equation ^ insures that the boundary conditions on 
the substrate (z = 0) are automatically satisfied for any 

Thus in what follows, we seek to solve for the 
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coefficients ^[^^ and 5^'^) using the equations that fol- 
low from the boundary conditions at the surface of each 
nanosphere, i.e. where rj = aj. Through rearrangement 
of the equations following from the boundary conditions 
on the sphere surfaces, one can eliminate the coefficients 
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. In Appendix A the linear set of equations deter- 
mining a\11 and B^^&ie derived [cf. Eq. (A4)]. 

In Sec. |III[ we present a series of numerical studies of 
plasmon resonance phenomena for nanosphere monomers 
and dimers placed on a substrate. To this end, we must 
solve Eq. ( |A4[). In order to do so, we truncate the summa- 
tions in Eq. (Isj) and also the equation system in Eq. ( A4 ) 
ail = L. The number of unknown coefficients in Eq. (KA) 
is then N{L -I- I)^ — 1. We use the same truncation 
limit for both the nanosphere-nanosphere interactions, 
as well as for the nanosphere-substrate interactions. The 
nanosphere-substrate interactions include both the inter- 
action of a given nanosphere with its own image, but 
also the images of the other nanospheres. Note that this 
is the first time, to the best of our knowledge, that all 
particle-substrate and particle-particle interactions have 
been taken consistently into account (to a given order). 
In several previous studies the interaction with the sub- 
strate has been taken into account to a high order, while 
the particle-particle interactions have been accounted for 
to dipolar or quadrupolar ovAe^^^^^. We shall see in 
the next section that the use of the dipole approximation 
(retention of only the terms with I — X) 'm the particle- 
substrate interaction is very inaccurate from a quantita- 
tive point of view save for the case when the nanospheres 
are quite far from the substrate. 

In passing we note that the formalism presented in this 
paper can be applied to extend the formalism used in 



In this section, we present a series of studies of the 
influence of a dielectric substrate on the plasmonic re- 
sponse of isolated nanospheres and nanosphere dimers. 
In addition, we find "hot-spots" created by plasmonic 
resonances between dielectric spheres and metallic sub- 
strates, as noted above. We shall also see that termina- 
tion of the hierarchy of equations at the dipole (L = I, see 
Ref. [121 or quadrupole {L = 2) order provides a very poor 
quantitative description of interactions between particles, 
and between the particles and the substrate. We remark 
that it is evident from earlier studies, which utilize a dif- 
ferent methodologjEHj that higher order harmonics must 
be included in the description of particle-particle inter- 
action, since the fields associated with "hot-spots" are 
highly localized around the points of nearest contact. 
Thus one must retain spherical harmonics to high order 
to describe these features. 

For the purpose of studying particle-substrate interac- 
tions, we first consider nanoparticles modeled by a di- 
electric function of the Drudc form^^ , 
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where wp is the plasma frequency and 7 is the inverse of 
the free carrier relaxation time. For the ambient material 
we have chosen vacuum, i.e. = I, and a dielectric sub- 
strate of £_ > 0. The virtue of model studies based on 
the form of Eq. ([5| is that we may choose the relaxation 
rate 7 sufficiently small so that much detail is evident 
in the calculated results. For the Drude model parame- 
ters we assume wp = 3 eV and 7 = 0.03 cV. After our 
discussion of nanospheres consisting of Drude metal, we 
present results for geometries incorporating silver (Ag) 
nanoparticles. Among metals which exhibit plasmonic 
response in the visible part of the optical spectrum, the 
damping rate in Ag is modest and numerous experiments 
employ Ag based structure^^. It should be remarked 
that the optical response of aluminium (Al) is described 
very well by the Drude model. Unfortunately the plasma 
frequency is very high, close to 15 eV, so the interesting 
plasmonic resonances in Al based materials lie well into 
the ultraviolet. In our view, it would be of great interest 
to see experimental probes of structures that incorpo- 
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rate Al nanoparticles, with attention to the appropriate 
spectral range. 

One possible indicator of "plasmonic activity" is the 
total dipole moment of one of our spherical objects. With 
p{uj) being the dipole moment of a nanosphere at angular 
frequency uj, we define the dimensionless dipole moment 
as 
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where a is the radius of the sphere in question, and Sq 
is the vacuum permeability. In terms of our expansion 
coefficients, Aim, the three Cartesian components of the 
dimensionless dipole moment are given by 
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Since the dipole moment in general is a complex vector 
quantity, the quantity we display in the figures below is 
the modulus of the total dipole moment given by 

p{u}) = \p{uj)\ = VF^, 

where f symbolizes the Hermitian transpose. In our stud- 
ies of the interaction of a single sphere with the substrate, 
we shall display the total dimensionless dipole moment, 
along with field enhancement factors for applied fields 
perpendicular to the substrate {Eq \\ z) as well as par- 
allel to the substrate (£^o II Moreover, for the dimer 
illustrated in Fig. [T] we shall present results for all three 
Cartesian components of the applied field. 



A. The Drude Monomer and Dimer 

We begin by considering a single Drude sphere in vac- 
uum (£_|_ = 1) located a distance h — 0.05a above a 
substrate. In Fig. [2] we present numerical calculations 
of the dimensionless dipole moment for three choices of 
the dielectric function of the substrate, £_ = 1,2, and 
10. The choice £_ = 1 corresponds physically to the case 
where no substrate is present. We present results for two 
choices of the applied field: (i) perpendicular to the sub- 
strate {z direction, red dotted curves) and (ii) parallel to 
the substrate {x direction, blue solid curves). Since the 
system is invariant with respect to rotation about the z 
axis, the response to an applied field parallel to the y axis 
is identical to that shown for x polarization. 

In Fig. [2j^a) , we present the response of the sphere in 
free space, from which the Mie resonance at the frequency 
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FIG. 2. (Color online) The dimensionless dipole moment, 
p(cj), for a Drude metal particle on a substrate of dielectric 
function (a) e_ = 1, (b) £_ = 2, and (c) £_ = 10. For the 
vacuum case, i.e. e_ = 1, we obtain the Mie result at huj = 
Tvjjpl^fZ ~ 1.73 eV. For all plots, we have h — 0.05a, up = 3 
eV, 7 = 0.03 eV, and L = 50. 



hujpj\/i w 1.73 eV is readily observed. While an isolated 
Drude metal sphere has a spectrum of multipole modes 

1/2 

at angular frequencies w/ — lop (Iji^ll + 1)) with I — 
1, 2, 3, . . ., only the dipole mode with Z = 1 is excited by 
an applied field whose wavelength is large compared to 
the radius of the sphere. In the presence of a substrate, 
higher order modes may be excited by a spatially uniform 
applied field, as we shall see. These will appear at higher 
frequencies than the Mie resonance, as suggested by the 
fact that > wi for / > 1. 

When the dielectric function of the substrate is £_ > 
£+, the spectral response of the nanosphere is altered 
significantly. First for the case of modest dielectric func- 
tion, £_ = 2, the (former) Mie resonance remains the 
dominant spectral feature and it is red shifted by the 
proximity of the dielectric substrate to the nanosphere 
[Fig. [2jb)]. A substantial splitting of the modes is ob- 
served when the response to a field parallel to the sub- 
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FIG. 3. (Color online) The dimensionless dipole moment, 
p(cj), for one of the particles in a Drude metal dimer, on top 
of a substrate of dielectric function (a) e_ — 1, (b) e_ = 2, 
and (c) £_ = 10. For all plots, we have h — 0.05a, d = 
0.1a, tjp = 3 eV, 7 = 0.03 eV, and L = 50. 



FIG. 4. (Color online) The plots show p{lo) for one of the 
particles in a Drude metal dimer for different values of h. For 
all plots, we have e_ = 10, d = 0.1a, ojp = 3 eV, 7 = 0.03 eV 
and L = 30. 



strate is compared to the response to a field perpendic- 
ular to it. We also see activation of a higher frequency 
mode. One might be tempted to associate this with ex- 
citation of the quadrupolar mode with Z = 2, but when 
the sphere is so close to the substrate, classification of 
the mode by the angular momentum quantum number is 
no longer accurate since a large number of / modes are 
mixed together. We require L, the cutoff used in the hi- 
erarchy of equations displayed in Eq. ( A4 1 , to be on the 



order of 30 (or more) to obtain converged results. 

If the substrate has a large dielectric function (e_ — 
10), then the response of the sphere is modified dramati- 
cally relative to the free-space case [Fig. [2]jc)] . The split- 
ting of the low frequency resonances, for parallel and per- 
pendicular excitation, is now very large. For both orien- 
tations of the applied field, the oscillator strength of the 
next highest mode is comparable to the low frequency 
(dipole) mode. We also see a third mode in the spectrum, 
so the symmetry breaking provided by the substrate now 



asserts itself prominently in the response of the sphere. 

The appearance of these higher order modes can be 
intuitively understood as follows. When we apply an 
electric field Eq to a nanosphere, it will generate local 
evanescent fields. When h is small, some of the evanes- 
cent fields are reflected from the substrate, resulting in 
a non-uniform field around the sphere. This causes the 
simultaneous excitation of many different I modes, mean- 
ing that the notion of discussing modes in multipolar ter- 
minology breaks down badly. The crosstalk between dif- 
ferent I modes is also the reason why we see higher order 
modes (e.g. "quadrupole" modes) in the dipole moment 
(p(a;)) of the spheres. 

We now turn our attention to a discussion of the re- 
sponse of a Drude dimer, as shown in Fig.[T] The radiuses 
of the two spheres are both assumed to be equal to a. The 
distance between the spheres is d = 0.1a, and they are 
both placed a distance h — 0.05a above the substrate. In 
Fig. [3] we depict the dimensionless dipole moment of one 
of the spheres in a dimer whose axis is parallel to the x 
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axis and, hence, to the substrate. 

When the dimer is placed in free space [Fig. [sf^a)] and 
the appUed field is perpendicular to the dimer axis, there 
is one dominant resonance. This is the Mie resonance of 
the single sphere, slightly blueshifted due to the particle- 
particle interactions. In addition, a second weak mode 
shows up at higher frequency. In contrast, when the 
dimer is excited by a field parallel to the dimer axis (blue 
solid curve), we see a sequence of collective modes red- 
shifted by large amounts from the isolated sphere Mie 
resonance. These results are in agreement with previous 
work on nanoparticle dimers in free spac^P^. As for the 
case of the single sphere on a substrate, the fields gen- 
erated from one sphere causes higher order modes to be 
excited in the other sphere, and vice versa. 

For a substrate with modest dielectric function [e_ = 
2, Fig.jsjjb)] we see a splitting between the dominant col- 
lective modes excited by a field parallel to z (red dotted 
curve) and that excited by a field parallel to y (green 
dashed curve). This is to be expected since the pres- 
ence of the substrate will break the rotational symmetry 
around the x axis. If the dielectric function of the sub- 
strate is substantial [£_ = 10, Fig. [3](c)] we observe dra- 
matic differences between the spectral response for the 
three directions of the applied field. The shift of the low- 
est frequency mode for II from the Mie resonance 
of the isolated sphere is particularly dramatic. 

These results demonstrate that placing metallic 
nanosphere dimers over a substrate with large dielec- 
tric function will give rise to substantial field enhance- 
ments. Also, strong dipole moment enhancements can 
be achieved over a spectral range very large compared to 
that realized for a single isolated nanosphere. Thus, as 
this example illustrates, the interaction between struc- 
tured nanoparticle arrays and a substrate of substantial 
dielectric function can allow one to design objects with a 
broad plasmonic spectral response. 

Figure |4] shows how the response of the dimer depends 
on the distance h above the substrate. A substrate di- 
electric function e_ = 10 was assumed in order to em- 
phasize the influence of the substrate on the response of 
the dimer. In Fig. |4|^a), where h = 2a, the spectral re- 
sponse is very close to that of the isolated dimer, shown 
m Fig. |3];a). We see clear interaction effects with the 
substrate when h = 0.3a [Fig. Qb)], but it remains true 
that the spectrum is qualitatively similar to that of the 
free dimer. The dimer has to be close to the substrate 
for the interaction effects to modify the spectrum even 
for the large substrate dielectric function used in these 
calculations [Fig.[4jc)]. 

We pause for a moment to comment on issues of con- 
vergence. It is common to employ the dipole approxima- 
tion to describe intersphere interactions, and also inter- 
actions of nanoparticles with substrates using the image 
method^. In Fig. [s] we present how the position of the 
lowest energy collective mode for the case Eq \\ z de- 
pends on L, which determines the number of unknown 
Aim coefficients in Eq. (A4). The frequently used dipole 




FIG. 5. (Color online) Position of the lowest energy resonance 
as a function of h, in the case where Eq \\ z. The blue dotted 
curve shows the result when the calculation is done in the 
dipole approximation (i.e. L = 1), while the red solid curve 
1/ = 30 show the conversed results. For all Drude 
dimer with d = 0.1a was assumed, and the substrate dielectric 
function was e_ — 10. 
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FIG. 6. (Color online) The square of the electric field, 
|B|^/|i/o|'^, at point 1 (blue solid curve) and point 2 (green 
dotted curve) as a function of frequency, for the Drude dimer. 
The results are for the case where Eo \\ x, and the substrate 
dielectric function is £_ = 10. The other parameters are the 
same as in Fig. [3] 



approximation corresponds to L = 1, and from Fig. [5] 
one observes that it is inaccurate even when the dimer is 
far above the substrate {h — 2a), and becomes gradually 
worse as h is decreased. The cutoff L must be on the 
order of 30 to obtain converged results for the parameter 
ranges explored in this paper. One may appreciate the 
reason for this from earlier worlP^l. When two spheres 
are quite close to each other, one encounters collective 
modes wherein the fields are concentrated in a small an- 
gular range near the points of closest contact. Similarly, 
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FIG. 7. (Color online) The field enhancement \E\/\Eu\ in 
the xz plane. At hu) — 1.35 eV the highest field enhancement 
is found between the sphere and the substrate; at hw = 1.45 
eV the highest field enhancement is found on the line con- 
necting the two spheres. As shown in the figure, Eq \\ x. The 
system parameters are £- = 10, h — 0.05a, d — 0.1a, and 
L = 50. Of particular interest is the area between the two 
lowest-frequency peaks, where the location of maximum field 
enhancement "flips" between the two points indicated. 



when one or more nanospheres are very close to a dielec- 
tric substrate, one encounters collective modes localized 
around the "south pole" of the spheres — the points clos- 
est to the substrate. One requires large values of the 
cutoff L if one wishes to describe such modes accurately. 
Notice, by the way, that the mode frequency is signif- 
icantly redshifted when the dimer comes very close to 
touching the substrate. In passing, we note that conver- 
gent results do not guarantee correctness of the calcu- 
lated potentials. In order to do so, one has to explicitly 
make sure that the boundary conditions are satisfied to 
the required accuracy at all points on all interface^^. 

In Fig.[6l we examine the nature of the enhanced fields 
in the Drude dimer at the two points indicated in the in- 
set. Again we have assumed d = 0.1a and h — 0.05a 
in these calculations. For most of the spectrum, the 



largest field enhancement is found at point 2, the hot- 
spot where the two spheres nearly touch. Notice, how- 
ever, that we have very large field enhancements also 
between the south pole of the spheres and the dielectric 
substrate, in particular at 1.35 eV [see Figs. |6]and[7]ja)]. 
In our view, the region where the dielectric substrate is 
very close to the bottom of the sphere acts like an effec- 
tive potential well which "traps" surface plasmons at the 
south pole. The surface plasmons sense the presence of 
the dielectric through the fields associated with them in 
the region outside the sphere. 

We can see from Fig. |6] that point 1 and point 2 
are "hot" simultaneously at roughly the same frequency. 
However, as one scans through a given resonance peak, 
near hoj = 1.4 eV in Fig. [6) the hot-spot moves from point 
1 to point 2 and conversely, depending on the precise 
value of the frequency. Thus we have another example 
of the phenomenon of the "moving hot-spots" discussed 
in a recent publicatiorP^. We illustrate this behavior in 
Fig. [Tj where we plot jiJl/lE'ol in a contour maiP3. A 
small energy shift of 0.1 eV is enough to change the shape 
of the field enhancement considerably, and move the hot 
spot from point 1 to point 2. 



B. Ag Monomers and Dimers 

The Drude model discussed in Sec. IIII Al is useful to 
examine, since one may model metals in which the plas- 
mons are damped very lightly. Thus one can explore 
detailed structure in the response of the model system. 
In practice, however, interest resides in realistic metals 
that display plasmonic response in the visible. In this 
respect, silver (Ag) and gold (Au) are the two metals 
most studied experimentally. While Au is indeed "plas- 
mon active" , the plasmons in this material are in fact 
rather heavily damped. Ag is a much better material in 
principle, even though in experiments oxide can form on 
its surface. 

This section is devoted to studies of the plasmon reso- 
nance properties of Ag monomers and dimers. Figure [8] 
shows calculations of the reduced dipole moment for a 
single Ag nanosphere placed a distance h = 0.05a over 
a dielectric substrate. For a free-standing Ag sphere in 
vacuum, the Mie resonance at fkv = 3.5 eV is readily ob- 
served [Fig. |8][a)] . The response of the sphere is modest 
fore_ = 2 [Fig.|8|;b)]. However, whene_ = 10 [Fig.jsjc)], 
we see a substantial splitting of the main resonance and 
activation of higher frequency modes occur. 

In Fig. |9] we show the response of an Ag dimer with 
its axis parallel to the substrate. As previously, we have 
assumed h = 0.05a and d = 0.1a. The strong interaction 
between the two spheres of a free-standing dimer in vac- 
uum can be seen from Fig.[9|^a) by noting the pronounced 
difference in response to an applied field parallel {Eq || x) 
or perpendicular {Eq \\ y or Eq \\ z) to the dimer axis. 
At least from the perspective of the dipole moment of 
each sphere, the influence of the substrate is not signifi- 
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FIG. 8. (Color online) The dimensionless dipole moment for 
an Ag sphere placed a distance h — 0.05o above a dielectric 
substrate of (a) e_ = 1, (b) e_ = 2, and (c) e_ = 10. The 
equation system was truncated at L — 50. 



FIG. 9. (Color online) The response for an Ag dimer close to a 
substrate, as described by the dimensionless dipole moment of 
one of the spheres. As in previous cases, d = 0.1a, h = 0.05a, 
and L = 50 were used. 



cant for S- — 2 [Fig.|9jb)], but we see substantial effects 
for the larger dielectric function, e_ — 10 [Fig. [9]Jc)]. 

While the dipole moment of the Ag spheres shows sub- 
strate effects to be weaker than those of the correspond- 
ing Drude monomer and dimer, the field enhancement 
effects are still substantial. When the spheres are ei- 
ther close to each other and/or close to the substrate, 
the resonances are highly localized in space and form 
so-called "hot-spots". This is illustrated by Fig. [TOj 
which shows the enhancement in the electric field inten- 
sity (|£Jp/|£^oP) for an Ag dimer. Hence, one can have 
local regions where the fields are strongly enhanced while 
their effect on the total dipole moment of the sphere is 
more modest. 

Regarding the field enhancement in the Ag dimer, de- 
picted in Fig. [TOj we see considerable enhancement be- 
tween the sphere and the substrate. This enhancement 
is caused by the proximity of the sphere to the dielectric 
substrate that creates a potential well where surface plas- 
mons can be trapped near the south pole of the sphere. 
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FIG. 10. (Color online) Intensity enhancement, |i?|'^/|iSo|'^, 
as a function of frequency of the applied field at point 1 and 
point 2, for the Ag dimer on a substrate of dielectric function 
e_ = 10. The remaining parameters are the same as in Fig. [9] 
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On resonance, the enhancement in the square of the field 
is close to 5 X 10^ [Fig. 10 . If one has SERS in mind, 



where the cross section is enhanced by roughly the fourth 
power of the field, then the Raman cross section would 
in this case be enhanced by 25 x 10^. Thus, the influence 
of the dielectric substrate on the enhanced fields realized 
for the dimer is very substantial. Although we observed 
full reversal of hot-spot positions for a dimer made from 
Drude metal, it appears that the larger attenuation of sil- 
ver (G[£j (uj)]) prohibits this phenomena in the Ag dimer. 
Hence, the dominant hot-spot is for all frequencies of the 
incident light located at point 2, in the gap between the 
two spheres. 
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FIG. 11. (Color online) The dimensionless dipole moment of 
a member of a dielectric dimer placed in close proximity to 
an Ag surface. The dielectric functions for the two spheres 
forming the dimer are both (a) ej — 2 and (b) Sj = 10. In 
both cases we have d — 0.1a, h — 0.05a, and L — 50. 



C. Dielectric Dimer on an Ag Substrate 

In the literature, primary attention is directed toward 
nano-scale objects fabricated from plasmon active met- 
als. We find that dielectric particles create localized 
evanescent fields that stimulate the formation of local- 
ized plasmons in a nearby metallic substrate. These sur- 
face plasmons are not excited on a flat metallic surface. 
In this section we consider a dielectric dimer with a fre- 
quency independent, real and positive dielectric function 
Ej placed close to an Ag surface. 

Figure [TT] depicts the dimensionless dipole moment of 
the dielectric dimers placed very close to an Ag substrate. 
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FIG. 12. (Color online) Intensity enhancements at two se- 
lected points for a dielectric dimer placed in the near vicinity 
of an Ag surface. The dielectric function of the spheres is 
e = 10, and other parameters are the same as in Fig. |11| 



As before, the separation between the two spheres is c? = 
0.1a and the height above the substrate is h — 0.05a. 
Dipole activity in the dielectric sphere is observed in the 
frequency range near the surface plasmon resonance of 
the Ag surface. Since the dielectric function of the sphere 
{sj) is frequency independent, this plasmonic activity has 
its origin in the Ag substrate. As expected, the effect is 
enhanced when the dielectric function of the spheres, Sj , 
is increased 



[Fig. 
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Figure 12 illustrates the frequency dependence of the 
intensity enhancement at a position between the dielec- 
tric spheres and the substrate (point 1) and between the 
spheres (point 2). Between the spheres (point 2) the plas- 
monic response of the substrate plays only a minor role in 
the intensity enhancement, whereas just below the south 
pole of the spheres, the plasmonic activity plays a more 
important role, producing higher intensity enhancement 
for a narrow part of the spectrum. The proximity of 
the dielectric spheres to the substrate converts the in- 
coming plane wave to an evanescent wave which excites 
surface plasmons in the Ag substrate. The consequence 
is that substantial intensity enhancements appear below 
the nanospheres near the surface plasmon frequency of 
the Ag surface. The physics is quite similar to the forma- 
tion of a hot-spot between a metal sphere and a dielectric 
substrate. 

The structure of the field enhancement near the plas- 
mon resonance is shown in detail in Fig |13| This fig- 
ure is qualitatively similar to as Fig. [7] Again, we are 
faced with the moving hot-spot phenomenon, since the 
hot-spot moves from between the spheres to below the 
spheres as the frequency of the incident light changes. 
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FIG. 13. (Color online) The field enhancement |i5|/|i5o| hi 
the xz plane for a dielectric dimer above an Ag substrate. At 
fiio = 3.02 eV the highest field enhancement is found between 
the spheres; at /Iw = 3.39 eV the highest field enhancement 
is found on the line connecting the two spheres. As shown in 
the figure, Eq \\ x. The parameters are Sj = 10, h = 0.05a, 
d — 0.1a, and L = 50. 



IV. CONCLUDING REMARKS 

We have formulated the theory of the interaction of 
non-periodic nanosphere arrays with a substrate, and for 
the case of monomers and dimers we have provided nu- 
merical studies of electric dipole moments induced by a 
uniform driving field, and also field enhancement gen- 
erated by excitation of plasmon resonances. While the 
focus usually is placed on the interaction of metallic 
nanoscale objects with metallic substrates or surround- 
ings, our emphasis has been on the interaction between 
metallic and dielectric materials: metallic nanoparticlcs 
on dielectric substrates, or dielectric nanoparticlcs on 
metallic substrates. Such systems are, in our view, better 
suited for experimental examination. 

For both these configurations we find "hot-spots" (i.e. 
local regions of high intensity) that are localized between 
the south pole of the nanosphere and the substrate. The 



physical origin of this behavior is an effective potential 
well created by the dielectric that traps and localizes plas- 
mons in the nearby metallic component. Consider, for 
example, a semi-infinite slab of a model metal described 
by the Drude model, Eq. ([5|. Let the metal lie in the 
half space z < 0, and let the half space z > be vacuum. 
The surface supports surface plasmons, and in electro- 
static theory these have frequency tup/y^ independent 
of wave vector. Suppose we instead fill the upper half 
space z > with a dielectric material whose dielectric 
constant is e+ > 1. This lowers the frequency of surface 
plasmons on the metal surface to wp(l 4- e+)~^/^. If we 
then imagine that the dielectric covers only a finite area 
on the metallic surface, clearly an attractive potential 
well is formed which can trap surface plasmons bound to 
the region where the dielectric is found. The frequency 
of these modes lies below the frequency band associated 
with those on the metal/vacuum interface. In our stud- 
ies, we have a rather different geometry. For instance, 
in one configuration we explored a dielectric sphere that 
is placed just a bit above the metallic substrate. The 
surface plasmons on the metal surface "sense" the pres- 
ence of the dielectric through their evanescent field that 
extends above the metal surface. We then find plasmon 
modes localized in the near vicinity of the south pole of 
the sphere. In the case of dimers, of interest is the "mov- 
ing hot-spot" phenomenon illustrated in Fig. [7] Small 
changes in excitation frequency results in a hot-spot that 
moves from one point in the structure to another. An 
earlier discussion provided an example of this behavior 
in a rather different structur^^. 

The hot-spots localized between metallic spheres and 
dielectric substrates, and also between dielectric spheres 
and metallic substrates, suggests that strongly enhanced 
non-linear optical studies may be carried out on diverse 
systems, not just those where all constituents are plas- 
mon active metals. It would be of interest to explore field 
enhancements not just for spheres placed near fiat sub- 
strates, but for other nanoscale objects of diverse shape 
as well. It should be possible to engineer structures in 
which large field enhancements are realized that can be 
exploited to study, for instance, adsorbates on insulating 
surfaces. 
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where Rij is the vector between the center of the sphere 
i and j, and 9ij and (f>ij are the polar and azimuthal 
angles, respecti vely , which describe the direction of Rij. 
In writing Eq. ( A2 1 , we have used 



Appendix A: Construction of the equation system 



The expression for the electrostatic potential, ipir), is 
shown in Eqs. ([T]) and ([s]). This expansion is of no use 

before one has determined the expansion coefficients A 



(i) 

and Bl'^. To do so, we combine the series expansion of 
and require the fulfillment of the boundary conditions 
at the surface of the spheres, i.e. the continuity of ip and 
ednip over any interface (where 9„ = n • V denotes the 
normal derivative ^I^H, Note that the boundary conditions 
at the interface at z = are already fulfilled through 
Eq. Q. For instance we may consider continuity of the 
electrostatic potential at the surface of sphere j. This 
gives the condition 



H{l,,mj,\k,m,) = V4^(-l)'"+'"^- 



21, + 1 



1/2 



(2/j- + l)(2/ + l) 

I + m\ / I — m 
li + rriij \lj + rrij 



(A3) 



1/2 



lim ipjivj) — lim < — r ■ Eq 



N 



N 



(Al) 



i=l 



1=1 



where / 
tion (t) 



and m = m 
denotes the binomial coefficient 



nij. Moreover 



the nota- 
The expansion 



described by Eqs. (A2) and (A3 1 can also be applied to 



i;i ± , 1 ■ ■ c 4- ■ the image multipoles located in the substrate (z < 0). 

where the notation a"^ means flj ±?7, where 77 IS mnnetesi- b f \ j 



mally small and positive. For this condition to be useful, 
and similarly for the equation following from the conti- 
nuity of the normal components of the displacement field 
D = — eeoV'0, we need to express the potentials ipi and 
^pi for i ^ j in terms of the coordinate system Sj cen- 
tered on sphere j. One can do this by using an identity 
employed by Bedeaux and Vliegei^^. This reads 



) 



(A2) 



^3 



We now have the electrostatic potential on each side 
of the surface of sphere j expressed in terms of the co- 
ordinates of system Sj. One may generate a system of 
equations for the unknown amplitudes by equating the 
coefficients of {9j,<f)j). When Eq. (Al) is combined 



with the condition that the normal components of the 
electric displacement field D should be continuous across 
the surface of sphere j, it is possible to eliminate the coef- 
ficients Bi'^ and to generate a linear system of equations 

that involves only When this is done, the following 
linear system of equations results: 



-&im,^: 



i,i,- 



i^j li.rrii 



R. 



Ij+U + i 



(A4) 



where /3 — (e+ — £_)/(e+ + £_) and Ij — 1, 2, 3, . . . , i expansion coefficients of the applied field £^0 in terms of 
and TOj = 0, ±1, ±2, . . . , ±Zj. The coefficients are the the spherical harmonics. The non-zero bim coefficients 
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(in the case of uniform Eq) are given h^^^^^ 



-Eo\l — cos Oq, 



bi±i = ±Eo\l — smt/oe 



Ti<Po 



(A5a) 
(A5b) 



where is the angle between the external field and the 
positive z axis and 4>o is the azimuthal angle which de- 
scribes the angle between the projection of the external 
field onto the xy plane and the positive x axis. As the 
Laplace equation is linear, we only need to solve for 3 dif- 
ferent directions of Eq {Eq parallel to x, y, and z). The 
response to an applied field pointing in any other direc- 
tion can be constructed through superposition of these 3 
cases. 



Equation ( A4) gives us N{L + 1) — 1 linear equations 



in the expansion coefficients A^^, and we have A^(L-l-l)'^ 
unknowns. The final equation results from the continuity 
of the normal component of D at the spherical interfaces. 
Taking the normal (i.e. radial) derivative of the Bqq term, 
we see that this term vanishes (9^.^ .Boo^o'(^i' ~ '^)- 
This means that Aqq = 0, related to the fact that the 
nanoparticles are assumed to carry no charge. Hence, 
the equation system is closed, and we can expect to find 
a unique solution. 

Finally, the resonances for an isolated sphere in a ho- 



mogeneous background of dielectric function can be 
obtained from Eq. (A4). By neglecting all contributions 

to keep 
(lAl, 



from other particles and image multipoles, i.e. 
only the first term on the right hand side of Eq. 
one is essentially left with the isolated sphere case. Un- 
der this assumption, the resulting equation can readily 
be solved to give 



Aim OC 



1 



le + e+{l + l)' 



where e is the dielectric function of the sphere. Hence, 
the resonance positions are determined by the zeroes of 
the real part of the denominator of Aim- 

^[le + e+{l + I)] =0. 

If we assume for e the Drude model with j — (e — 
1 — Wp/w^), which in our case is a good approximation, 
we get the following resonance frequencies for the isolated 
sphere: 



UJl = UJp 



I 



2/-M 



For systems containing more than a single, isolated 
sphere, such as the ones discussed in this paper, these res- 
onance frequencies are typically modified due to particle- 
particle or particle-substrate interactions. 
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